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METHODS OF SEPARATION OF VARIABLES IN TURBULENCE THEORY 


by Shunichi Tsuge 

Nielsen Engineering & Research, Inc. 


SUMMARY 


Two schemes of closing turbulent moment equations are proposed both 
of which make double correlation equations separated into single-point 
equations. The first is based on neglected triple correlations. Grid- 
produced turbulence is described in this light as time-independent, 
cylindrically-isotropic turbulence, showing the existence of two modes 
in fluctuations, one obeying and the other disobeying Taylor's hypothesis 
respectively. Application to wall turbulence guided by a new asymptotic 
method for the Orr-Sommerf eld equation suggests existence of a neutrally 
stable mode of essentially three dimensional nature. The second closure 
scheme is based on an assumption of identity of the separated variables 
through which triple and quadruple correlations are formed. The result- 
ing equation adds, to its equivalent of the first scheme, an integral of 
nonlinear convolution describing a direct energy-cascading. 



I . 


INTRODUCTION 


A. Separability of Double Correlation Equations 

The first sound basis on which equations governing double and 
higher-order turbulent correlations are formulated is due to Karman and 
Howarth (ref. 1). Among other correlation or Reynolds' stress equations 
constructed in various ways on the common basis of Navier-S tokes frame- 
work this equation has an additional kinematical attribute that the 
Karman-Howarth formalism is consistent also with the BBGKY hierarchy 
(ref. 2). In other words, it alone is concurrent with Liouville’s 
equation; a master equation of nonequilibrium statistical mechanics. 

This fact has been revealed in reference 2 by starting with the two-point 
BBGKY equation transformed into fluid-dynamic moment equations through 
13-moment expansion, leading for incompressible turbulence to an equation 
identified a posteriori with 

< A a • (NS ) 0 + Aa. (NS ) . > = 0 (1) 

3 l l 3 


where we have defined the operator (NS)^ ky 


(NS) 


+ u 

dt 


r dx 


1 ^P 

ax 2jU ' p 3x. 


u . + 


( 2 ) 


and an averaging procedure by the symbol < >. In these expressions 
quantities without and with caret (~) denote those at space points x and 
x, respectively, cm and p are the fluctuating fluid-dynamic velocity and 
pressure, respectively, and Am is the instantaneous velocity fluctuation 

Aa. = a. - u, (3) 

3 3 3 


namely, instantaneous deviation around the average velocity 

u . = < u . > ( 4 ) 

3 3 

Equation (1) gives us a 'recipe' of how to cook fluid-dynamic equa- 
tions correctly for the correlation equation to be consistent with the 
BBGKY formalism. It requires that the velocity fluctuation to be 
multiplied on the Navier-S tokes operator is at a different space-point 
and at the same time, so that it prohibits using the same space-points 

x = x in Eq ■ (1) in contrast to what is often done in constructing model 

~X 

equations. Because of the seven-dimensional structure in (x,x,t) space 


2 



this equation, known to f luid-dynamicists without its having the firmer 
basis recognized, has been thought as "totally intractable" (ref. 3) for 
flows with more realistic geometry. Applications, therefore, have been 
limited so far to homogeneous turbulence, where physical quantities 
depend only on relative location of the two points. If, however, an 
assumption is made that triple correlation in Eq . (1) is negligibly small, 

the equation is seen to be separable independently of its flow geometry. 
Actually this observation has rendered the equation separated out of (1) 
an Orr-Sommerf eld equation (ref. 2). 

The feasibility of the variable separation has been found originally 
by a heuristic observation, therefore, it has been believed for some time 
that such a separability would be restricted to incompressible flows. 
Actually, allowance for the density fluctuation in the case of compress- 
ible turbulence makes it difficult to judge its feasibility on intuitive 
basis. We will see, however, that the separability of the equation is 
examined more easily through the kinetic equation rather than its moment 
(fluid-dynamic) versions. Following this guideline we will examine the 
structure of the kinetic equation of the second BBGKY hierarcy governing 
a correlation function between two phase-space points z = (x,v) and 
z = (x,v) . The correlation function as meant here is defined by 


g(z,2) = < A A$> (5) 

where denotes instantaneous fluctuation 

A< M - f (6) 

of the microscopic density (ref. 4) around its average 


f = <i> 


(7) 


namely, the Boltzmann function. There are direct relationships between 
the kinetic mean and fluctuation variables given here and those in fluid- 
dynamic space as defined loosely by (3) and (4). These are given con- 
sistently with the classical kinetic theory as 



6 dv 
A ^ dv 


( 8 ) 

(9) 


where p is the density, and where the ad hoc assumption of incompress- 
ibility is made for brevity of the expessions. By taking product of the 
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instantaneous velocity fluctuations (9) at different space points x and 
x, and then taking average, we have an expression for turbulent velocity 
correlation tensor 


<AUjAu^> = 


-2 


v . v p g (z , z) dv dv 


( 10 ) 


where the definition (5) of kinetic correlation function g has been 
incorporated. The function g of form (5) includes, by its nature, two 
types of fluctuations in gases; namely, the thermal and the turbulent 
fluctuations (ref. 2). In what follows, the first effect is ignored 
since it is by far the smaller in magnitude, also we are interested only 
in fluid mechanics. If g = 0 the gas is in (binary) molecular chaos. 


B. A Perspective of the Work 

Examination of the structure of the two-point kinetic equation 
reveals a fact that the variables are separable into two groups (z,t) 
and (z,t) under any conditions of geometry or compressibility of a flow 
if, in the equations, the term of triple correlation is negligible com- 
pared with those of binary correlations, 

<Atf A l A<$> = 0 (11) 

in other words, if the gas is in " ternary" molecular chaos. (First 
separability condition, section II.). The separated equation in (z,t) 
space forms a closed set together with the Boltzmann equation in which 
g of (5) is retained in the collision integral. Each of the coupled 
equations is moment-expanded by means of 13-moment method, providing 
with 26-moment fluid-dynamic equations. Coupling between mean and 
fluctuating quantities is effected through only two terms representing 
the Reynolds stress and the turbulent heat flux density in the mean 
momentum and the energy equations (section III.). This method gives us 
a means for a more direct formulation of grid-produced turbulence as 
time-independent and cylindrically-isotropic phenomena, enabling us to 
find a solution hidden behind the conventional assumptions. Also the 
method provides us with dispersion relationships, showing that the grid 
turbulence consists of two modes, the one stationary relative to the 
grids, the other frozen to the fluid; the latter alone obeys Taylor's 
hypothesis. An initial-value problem is posed in such a form that com- 
plete knowledge of the velocity correlation is provided in terms of 
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frequency-analyzed transverse correlation data prescribed at an initial 
plane (section IV.). Application of the proposed formalism to incom- 
pressible wall turbulence is also undertaken. The separated fluctuation 
variable is, in this case, governed by Orr-Sommerf eld r s equation (ref. 2). 
Difficulties in solving the equation, caused by anomalous velocity 
variation near the wall is circumvented by devising a new asymptotic 
method. Essential role of three-dimensionality of fluctuations is main- 
taining Reynolds stress is elucidated using the tool (section V.). In 
section VI. the first closure condition is discarded and is replaced 
with the second closure condition 

+ + + 

<p = * TT (12) 

"f* ^ 

where (p and <p denote, respectively, separated variables through which 
triple and quadruple correlations are formed. This closure retains terms 
of the triple correlations in the two-point equation, deducing this and 
three-point equations in an identical form through the variable-separa- 
tion. The resulting separated equation describes contribution from the 
triple correlation as a convolution integral, indicating primary role 
of the term to be energy-cascading that is far reaching in the frequency 
space . 
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SYMBOLS 


A 


o 

B ( f ) 


B 

o 

b 


i (i?l,2,3,4) 
f 

rf(z> 


f T1 (z,z) 


g (z , z) 

H. (n) 

ID 


i (i = 2 f 3) 
J (z | z) 
T (n) 


k 

L 


(0) 

( 1 ) 


magnitude of total wave number defined by ( 101 ) 
complex amplitude of Laplace mode defined by (82) 
complex amplitude of Oseen mode defined by (82) 
quantity defined by (B.13) 
operator defined by (124) 

complex amplitude of Oseen mode defined by (81) 

velocity gradient at the wall nondimensionalized by 
uniform flow velocity u and geometrical characteristic 
length l 

isothermal speed of sound (39) 
wave velocity defined by ( 101 ) 
differential operator d/dy 

four independent solutions of Orr-Sommer f eld equation 
given by (115) 

Boltzmann function 

microscopic density, or instantaneous number density in 
the phase space 

defined by ( 6 ) 

two-particle distribution function 
quantity defining $2 through (B.15) 
correlation in the phase space 
threedimensional Hermite polynomials 
integrals given by (B.30) 

Boltzmann collision operator defined by (15) 
quantity given by (B.36) 

complex wave number in stream direction 
= (k^,k 2 ^ 3 ) wave number vector 
operator defined by (54) 
linear operator defined by (113) 
linear operator defined by (114) 
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SYMBOLS (Continued) 


m 

m . 

3 

(NS) j 
P 
P 
P 


jk 


Q ■ 


q j 

„ (n) 

q j t. . . 

R 

R 

r (J ' K) 
i j . . ,£m . . 

r 

S 

u 


u . 

3 

u . 

3 

A u . 
3 


III 


IV 


x 

Y 


geometrical characteristic length 
mass of a particle 
mean mass flux density 
nonlinear operator defined by (2) 

mean pressure; defined by (40) for compressible flows 
instantaneous pressure 

turbulent stress deviator tensor defined by (46) 

turbulent heat flux density defined by (47) 

separated velocity fluctuation variable for incompressible 
flows 

expansion coefficient of (J) defined by (48) 

Reynolds number defined by (79) 
designating real part 

turbulent correlation in physical space 

defined by (84) 

quantity defined by (B.16) 

velocity of a uniform flow 

j-component of mean velocity 

j -component of instantaneous velocity 

j -component of instantaneous veloicty fluctuation 
defined by (3) 

quantity defined by (B.7) 

coordinate in the (molecular) velocity space 
quantity defined by (B.7) 

Wronskian given by (B.23) 

Wronskian given by (B.27) 

Wronskian given by (B.34) 
coordinate in the physical space 
auxiliary variable defining Y ^ 
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SYMBOLS (Continued) 


Y . 

J 

Y (n) 

Z 

z 

a 

r i ( i ; 2 , 3) 

Y 
6 
£ 

(i;l,2,3,4> 
4> (z) 

n 

A 

X 

M 

Y 

P 

T 

0 

X 

tr 

n 

n 

a 

(jO 

u)(5) 


defined by (95) and (97) 

n-th order eigenfunction in the successive approximation 
(109) obeying (110) through (112) 

auxiliary variable defining Y. through (97) and obeying 
( 100 ) 3 

z = (x,v) 

wave number in streamwise (x) direction (defined by (95)) 
quantities defining through (B.21) 

wave number in spanwise (z) direction (defined by (95)) 

designating deterministic perturbation 

scaling parameter defined by (105) 

four independent solutions of (107) 

separated fluctuation variables 

defined by (70) 

linear operator defined by (133) 

thermal conductivity coefficient (in the main text) . 
ah hoc variable defined by (B.5) (in Appendix B) 

viscosity coefficient (in the main text) 

ah hoc variable defined by (B.5) (in Appendix B) 

kinematic viscosity 

density 

quantity defined by (B.13) 

threedimensionality of the fluctuation defined by (98) 
defined by (24) 

nondimensional molecular velocity defined by (38) 
independent variable defined by (105) 
operator defined by (18) 
operator defined by (127) 

constant separating variables (see (16) and (17)) 
local Maxwellian distribution defined by (37) 
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SYMBOLS ( Concluded ) 


< > taking average 

~ symbol designating a second point 

symbol designating a third point 

Superscripts 

* complex conjugate 

+ pertaining to quantities nondimensionalized by the 

friction velocity and the kinematic viscosity 
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II. 


KINETIC EQUATION IN SEPARATED FORM 


A. BBGKY Hierarchy Closed at Two-Particle Equation 

According to the standard BBGKY hierarchy formalism the one- and 
two-particle equations with the latter being truncated by closure con- 
dition (11) of ternary molecular chaos, are written, respectively, as 


3 

at 


+ 



J(z|2)[f(z)f(z) + g(z,z)] x= £ * 


a 

at 


+ 


V 



-y 

V 



g (z, z) 


J(zjz) [f (z) g (z , z) 


( 13 ) 


+ f (z)g (z, z) ] x _~ + J(z|z)[f(z)g(z,z) 


+ f (z)g (z,z) ] - = ~ 


(14) 


where f is the Boltzmann function defined by (7), g is the kinetic 
correlation function given by (5) and connected with two-particle dis- 
tribution function f^ of the BBGKY formalism through 

g ( z , z ) = <Atf A^> 


= f IX (z,z) - ff 

In this formula and in what follows an abbreviated expression f = f(z) is 
employed. The operator J(z|z) in Eqs. (13) and (14) is the Boltzmann 
integral operator which acts on the field particle z in the manner 


with 


J(z j z) [g (z) ] x= £ 


i r 
m 

J 


[g (z 1 


z ' ) 


dK = 2rr I v - v | bdb 


~ if 

g (z , z) ] dKdv 






(15) 


In order for the two particles z and z to execute a molecular collision 

-y- ■+ 

we need to impose a condition x = x, which is indicated as a subscript 
to the collision integral. In this collision integral the quantities 
with prime denote the state before a collision leading to unprimed state 
upon collision, b is the impact parameter, and m is the mass of a 
particle . 
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B. Separation of Variables 

The form of equation (14) suggests possibility of the six-dimen- 
sional equation to be separated into two groups of variables z and z by 
assuming 


g(z,z,t) = <Ajj A|j> = 


<j)(z,t;{jj) <(>(z,t;u>) dto 


(16) 


This expression, substituted into (14) and rearranged properly leads to 
the equation of the following form 



4 ) 


$ 


ico 


(17) 


with 


<f>(z) - J(z|z) [f (z)<Mz) + f(z)<Mz)] x= ~ U8) 

and with ft defined likewise, namely, with variable z replaced by z. Note 
that the time t is the only variable common in both ft and ft. In line 
with the general separability theorem, the left-hand side of (17) is 
only a function of z, whereas the middle term depends only on z; there- 
fore ico is necessarily a constant separating the variables. The govern- 
ing equations for <J> and $ take, accordingly, the following forms 

ft$ - iuxf) = 0 
ft$ + icotj) = 0 

In view of expression (18) we see that and <$> obey a linearized Boltzmann 
equation in a generalized sense (co ^ 0) . It is also seen that the 
separation constant co has the dimension of the frequency. It would be 
natural to postulate that correlation function g be symmetric with 
respect to its arguments (z,z) . As is easily seen from (16), this con- 
dition is met either with $ = cj> or with 

$ = 4 >* ( 20 ) 

where the symbol (*) denotes the complex conjugate. This latter require- 
ment is satisfied, in view of Eqs . (19) and reality of operator ft, if 

and only if 

to* = co , or to; real (21) 



ft cj> 


3 

at 


+ v 
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In contrast, the former condition (0 = 0) together with Eqs . (19) 

requires necessarily w = 0, resulting in a trivial solution (ref. 2) 

0 ~ 0 . These observations indicate that 0 is essentially complex. 

It seems to be worthwhile to remark at this point certain parallelism 
of function 0 to Schrodinger ' s function V, and of Eq . (19) to Schrodinger ' s 

equation 


HY - ifidV/dt = 0 


( 22 ) 


where H and are Hamiltonian operator and Planck's constant, respec- 
tively. Both of the equations have the imaginary factors associated with 
temporal operators based on essentially physical basis, not as a result 
of mathematical method of solution (c.f. Eq. (35) below). Also, both of 
the dependent variables 0 and ¥ are related through the same bilinear 
forms to the same quantity, namely, the probability density. Its actual 
form in the case of the Schrodinger function is well-known: the quantity 

M / * l Fdx gives the probability of finding a quantum particle within a three- 
dimensional volume dx in the physical space [x,x+dx]. The counterpart 
expression for 0 follows from (16) and (20), namely, 


g(z,z) = <A^ Aj$> = 


X (z , z ; u>) dw 


(23) 


with 


x(z, 2 ? 0 )) = R{0 (z?u>) 0* (z;ai) } 


(24) 


where the symbol R{ } denotes taking the real part. In fact X is bound 
to be a real quantity because expression (23) defines its physical 

meaning as the spectral density (in the frequency-space) of turbulent 

a -y X 

correlations in the phase space. More precisely, expression <A^ A^>dvdv 
gives the probability of finding a particle located at x in a cell 
[v,v+dv] of the (molecular) velocity space and another one that is 
located at x, interacting with the first one, to be found in a cell 
[v,v+dv] . Therefore, the quantity 

R { 0 ( z ; oo) 0* (z ; a)) }dv dv dw (2 5) 

gives the probability of finding such two particles in the frequency 

-+ -£ 

subspace [oj , co+dco] , in addition to their being in the cells dv and dv, 
respectively. Since the function 0 is attached physical meaning only 
via (25) as in its equivalent H^tdx for the Schrodinger function, the 
following transformation 
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( 26 ) 


j ia^ 

<$> ^ e <J> 

where a stands for any real quantity, keeps observable quantities 
invariant as in quantum mechanics. 

Expression (25) constitutes a basis on which any two-point quanti- 
ties of fluid turbulence can be calculated. For example, the turbulent 
correlation between velocities at two points x and x is, from (23) and 
( 10 ) , 


<AUjAu^> = R{Jdio( 


v j 0dv) ( 


v^$*dv) } 


(27) 


and likewise for any other macroscopic quantities definable in the 
classical kinetic theory. 

With these physical backgrounds the Boltzmann function f and the 
separated fluctuation variables $ form a closed set consisting of Eqs . 
(13) and (19) . These are written, in view of (23) and (24) , as . 


f 

a 

at 


+ 


ax 


— f = J(z|z) [ff + R{ <Mz)<f>*(z)dw}] 


x=x 


-id) -f — t V 


_a_ 

a^ 


j(z|z) [f (z)<f>(z) + f (z)<Mz)l 


x=x 


(28) 

(29) 


Needless to say that Eq . (28) is exact in the present context, whereas 

Eq. (29) which is a result of separation of variables is subject to, 
but only to assumption (11) of ’ternary 1 molecular chaos . We see, there- 
fore, that the separability is not affected by compressibility or geom- 
etry of the flows, proving the point at issue stated at the beginning in 
Section I. A. 


C. Implications of Dual Time-Dependency of Eq. (29) 

As a result of introduction of the separation constant w having the 
dimension of the frequency, Eq. (29) is seen to have dual time dependen- 
cies through explicit (3/3t) and implicit (-ioj) temporal operators. We 
will show that the former acts on deterministic unsteadiness associated 
with its mean motion, whereas the latter reflects stochastic unsteadiness 
of the turbulent motions that is averaged out, so is not to appear 
explicitly, in the one-point regime (Eq. (28)). Regarding the explicit 
one an explanation will be given through an assertion that if f is 
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time-independent, so is <J>. The contrapositive of the statement is; if 4> 
turns out to be time-dependent, it is caused by time-dependency of f. In 
other words, explicit time-dependency (3/3t ^ 0) is due only to unsteadi- 
ness in mean motion of the turbulent flow. 

Proof of the assertion, in its first expression is the following; 

For a time-independent f, Eq . (29) allows a formally time-dependent 

solution 


4>(z,t) = exp (-ito Q t) <!> o (z) (30) 

The invariance, however, of physical quantities under transformation (26) 
applied here assures that the temporal factor in (30) does not affect on 
physical quantities, so does not generate unsteadiness of a new class to 
be superimposed on the steady background flow. (Note that of (30) be 
real in order to keep g(z,z) of Eq . (23) symmetric with respect to z and 

z. ) 

We have thus proved that for steady mean flows temporal factors 
are involved in $ only through implicit dependence (-iu>) . This is a 
fundamental property of turbulent fluctuations as distinct from 
deterministic perturbations such as Tollmien-Schlichting waves that 
are bound to be explicitly time-dependent. 

It seems to be difficult to draw a clearcut picture as to the 
implication of the implicit time-dependency at this stage. In a sub- 
sequent section (Section IV. C), however, the quantity co introduced in the 
theory as a separation constant will be given a firmer physical basis 
through a dispersion relationship derived for grid turbulence, where 
oo/k^ (k^; the wave number in stream direction of the turbulent fluctua- 
tion) is interpreted as the propagation velocity of a turbulent vortex. 

In this connection the implicit time-dependency may be thought as having 
its physical origin in stochastic elements of coherent turbulent vortices; 
for example, the lifetime of retaining their identity as observed in a 
free shear layer (ref. 5), or stochastically repeated up- and down-washes 
in the wall-shear layer considered as being caused by aligned longitudinal 
vortices (ref. 6). 

D. Turbulent Fluctuation <J) as Distinct from Deterministic 

Perturbation 5f 

Suppose a small time-dependent disturbance 6f(z,t) superimposed on a 
steady basic flow described by the Boltzmann function f which is time- 
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independent otherwise. The small disturbance 6f as meant here is of 
fluid-dynamic nature, regardless of its kinetic-theoretical expression, 
as would be seen from its definition 

13 

6f = l (3f/3A,)6A, (31) 

A=1 A A 

where A-^ stands for 13 mean fluid-dynamic variables specifying the 
Boltzmann function (ref. 7) . 

Let us observe turbulent phenomena in which the two fluctuations Sf 
and <p are coexistent. One such example is a boundary-layer undergoing 
natural (ref. 8) or forced (ref. 9) transitions where Sf corresponds to 
naturally or artificially generated sinusoidal perturbations. Under 
this condition the function (j) is necessarily responsible for turbulent 
part of the flow observed as bursts and/or turbulent spots in the experi- 
ments. In these experimental observations where both deterministic and 
turbulent fluctuations are operating, we are especially interested in 
the situation where condition 

o ( I 4> I ) ~ O ( | 6 f | ) 

holds. For superimposed perturbation onto a basic distribution which, 
otherwise, is time-independent; 

f -* f + Sf 

we see from Eq. (28) that f and <5f obey the following equations; 


3 

9t 


-> 

+ v 





f = J (z | z) [ff + 6f 6f + R{ 44*du)}] 
6f = J (z | z) [f 6f + f 6f ] x=5 


x=x 


(32) 


(33) 


where term 6$ induced responding to 6f is ignored because of its second 
order smallness. Since Eq . (33) is linear in 6f and depends on t only 

through the derivative it reduces, by putting 

6f = R{exp ( -iu> o t) 6h} (34) 


to 
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( 35 ) 


' - a 1 

-iu) + v — <Sh = J(zlz) [foh + f6h] 

l ° 3xJ x=x 

The requirement for reality of 6 f in (34) comes from the fact that the 
f luid-dynamic perturbation 

6A = a (z) 6fdv 

where A is a macroscopic variable and a is its particle-dynamic equivalent, 
is an observable quantity such as the Tollmien-Schlichting wave. 

It is noteworthy coincidence that Eqs. (29) and (35) governing <J> and 
6h, respectively, are identical in form under steady flow condition 
(9/9t = 0) in the former. The coincidence, however, breaks down at the 
point of their different ways of interaction with the mean flow, reflec- 
ting their difference in physical nature. It is readily seen from Eq . 

(32) which reads 

+ V ~ f = J(z|z) [ff + R{exp (-iw t) x6h}R {exp (-iw Q t) 6h} 

’■ 9t 3x' 

» 

+ R{ <H*dw}] ~ (36) 

j x-x 

where Eq. (34) has been used. An obvious distinction is generation, on 

one hand, of (temporal) higher harmonics, namely, of the factor 

exp(±2iw o t) to appear in the quadratic term of the deterministic perturba- 
tion 6h, whereas the temporal factors cancel out from the term $0* of the 
turbulent fluctuation as discussed in Section II. C. 

This feature of nongeneration of higher harmonics, spatially or 
temporally, has been originally found in the physical space for the 
Reynolds stress (ref. 10), in contrast to the similar terms of nonlinear 
stability theory (refs. 11,12) which give rise to these effects. It will 

be seen (Section III. A) that the Reynolds stress in the physical space 

emerges from the bilinear term of Eq . (36), proving the identity of the 

two assertions. A direct observation by Roshko (ref. 5) on coherent 
vortices in a free turbulent shear layer exhibits evidence of no genera- 
tion of such higher harmonics in spite of their characteristically non- 
linear patterns. 
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III. A CLOSED SET OF FLUID EQUATIONS 


A. Review of One-Particle Turbulence Formalism 

Fluid-dynamic version of the 'turbulent' Boltzmann equation (28) has 
been derived in ref. 13 using the 13-moment expansion carefully extended 
to turbulent cases: We expand the Boltzmann function in the following 

form 


with 


f(z) = p co(£) 


w(S) = 


1 ♦ £ "8 + 


Q • 


5pc j 


<3) 




(2 TTC 


2, 3/2 


exp 




(37) 


J 


where p is the average density, and H's are the Hermite polynomials 
defined in terms of nondimensional velocity variable f given by 


-1 

v . - p m . 

= -J 1 

3 C 


(38) 


In order to preserve the original formality of this scheme (ref. 7) , it 
is requirt 
chosen as 


is required that m_. be taken as the average mass-flux density and C be 


C = (p/p) 


1/2 


(39) 


where p is a 'pseudo' pressures, distinct from the average thermodynamic 
pressure n by the added contribution due to turbulence; 


P - * + (3p) _1 r! 1 ' 1 ’ 

D * J 

For exact definition of the term ^ see Eq. (41) below. 

J t 3 

On the same background the kinetic correlation function g which 

appears in the collision integral of Eq. (28) is expanded in a double 

~y ~ 

series of H(£) and H = H(£) as 


(40) 


g(z,z) = w(£;)a)U) 


(3,3) 

I 


r (J,k) 
i j * « « » 


( J , K) = ( 0 , 0) C J C K J!K! 


, (J) 
ij • 




(41) 


A word must be mentioned at this point of truncating expansions (37) and 
(41) at third-order Hermite polynomials in line with the 13-moment regime. 
In laminar (classical) cases this is perfectly legitimate on fluid-dynamics 
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point of view, because the neglected terms are of 0(e) (e: the mean free 

path), no information on fluid-dynamics are lost by this truncation. For 
turbulent cases, on the other hand, the neglected terms include fluid- 
dynamic quantities of 0 ^ 2 )# namely, terms of double correlations 

A 2 - <AA AB> 

where A and B are any gasdynamic variables. Neglecting these terms, 
however, is justified only for Maxwellian molecules because these terms 
are not to intervene in the gasdynamic equations due to a special 
structure of the collision integral peculiar to this type of molecules 
(ref. 6 ) . 


Once the moment expansion for f and g have been assumed in forms 
(37) and (41) fluid moment equations out of one-particle kinetic Eq . (28) 

are determined uniquely through a standard procedure, leading to 


~ 3 m. 

3p , k 

3 1 3x, 


= 0 


3m . 

1 

3 t 


3x. 


rim m . 

Lj. 


+ tt 5 + p .. , 
Dk *jkJ 


= 0 


(42) 

(43) 


with 


_3_ 
3 1 


1 


+ e 


3x, 


m i m k , ^ 


h m . 
+ -JL 


[ 2p‘ 


P P jk 


+ Q, 


= 0 


(44) 


e = | -rt + R^ 1,1 ^ 

2 2p r , r 


h n e + 7T , 


p = (p ) + I rU.D 

P jk P jk NS + p R j,k 


J k = (c *k } F + 6^ R k 


( 1 / 2 ) 


(45) 

(46) 

(47) 


where (p., and (q. ) ^ are the viscous stress tensor and the heat flow 
^jk NS k F 

vector connected with the average gasdynamic quantities through classical 

Navier-Stokes 1 and Fourier 1 s laws, respectively. Note that turbulence 

corrections to the stress tensor and the heat flow vector comprise only 

one term in each case; the Reynolds stress p and turbulent heat 

(12) 3 t K 

flux density (5/6 ) R^ ' . If the gas is not monatomic, the factor 5/6 

is to be replaced with Y/3(Y-1), where Y is the adiabatic index. A 
noteworthy remark at this point is that all the equations (42) through 
(47) are written in compressibility-invariant form, so that the transport 
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relationships (46) and (47) , each expressed in a form with single-term 
turbulence correction, is valid for compressible flows as well with no 
recourse to mass averaging (ref. 13). Also, deduction up to this point 
is exact within the framework of fluid description for Maxwell molecules, 
independent of any closure hypothesis to be effected at a subsequent 
hierarchy. We note that the key to the compressibility-invariant expres- 
sions of turbulence equations is to employ variables proportional to 
density, namely, the mass-flux density and the pressure to replace the 
velocity and the temperature, respectively. With this code in mind, all 
the equations derived in this section are reproducible using a priori 
gasdynamic equations (ref. 13). 


B. Fluid Moments from Separated Two-Particle Equation 

Our task in this paragraph is to derive from Eq. (29) equations 
governing the Reynolds stress and the turbulent heat-flux R^ 1,2 ^ 

appearing in turbulent transport relationships (46) and (47), respectively. 
For this purpose we expand the dependent variable <J> in the form: 


<J> (z) = 0 ) (O 


.(o) + h 


(i) 


Hf 1) + 


(2) 

2C 2 


h ( 2) + it 

3 k 5(r 


(3) 


(3) 


(48) 


This form of expansion assimilates the 13-moment expansion (37) for the 
Boltzmann function except that the coefficient of the first order Hermite 
polynomial is nonvanishing, also the tensor coefficient of the second 
order has a nonvanishing trace; 



(49) 


This situation arises from the fact that we have no degrees of freedom 
left out to suppress these quantities because these have been depleted 
in defining the expansion form of the one-particle distribution. 
Expansion coefficients q’s of (48) are related to those R’s of (41) 
through (23) and (24). For example, the Reynolds stress and the turbu- 
lent heat flux are given in bilinear forms of q's through putting x = x; 


r 11 ; 11 = r{ 
3/* 


qj 1 ) (x) q 


, (1) * (x) dco} 


r (1 ' 2) 

1 




q U) (X) 
D 


q ^ 2 ^ (x)du) 



50) 
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Equations governing q's in the fluid-dynamic space are obtained from 
(29) through the standard moment expansion procedure familiar in the 13- 
moment theory and are detailed elsehwere (ref. 14) . Only the final form 
of the equations would suffice for what follows: 


(1) 

^ <0) * TF- 


= 0 


(51) 


^ 3 


9x _ 


r 


9x 



r „ (i) 

3 

q 

9x . 

p 

3 

~ -V 

l ) 


2 . 3 

7 6 j r 


+ q 


(0) 


Dw . 


= 0 


(52) 


3 ( 2)1 


L hr p 


+ q. 


:i) 


* p (2 » 

o w 

r 

9 

3x 1 
r 

5 r 2 n 

2 C q 

Dw 

r 

9w . 

9 

i — 1 

— u\ 
t r 1 

Dt 

y 9X 1 
r 

9x . 
L 3 

l p J 


9 

» r^’fl 


3 

^ 3 fq ,2, ll 


R 3x 1 p 

- r 

\ ✓ 


rq 


3x 


( 1 ) 

V) 


- ^ 6 , 


3 jr 3x, 


= 0 


(53) 


with the following definitions of the symbols, 

Lq i - iuq + || + „ r q 


W = iq (2 » + C 2 q (0 > 


m 

w = 
r P 


Dt 9t w r 9x 


(54) 

(55) 


where u and X denote the molecular viscosity and the molecular heat con- 
ductivity , respectively . 
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The system of Eqs. (51) through (53) comprises five equations to 
determine five unknowns q (0) , q^ 1 * , and q (2) (or p (2) ) that correspond 
to <5p , pSw . and pR6T (or <5p) , respectively, of deterministic perturba- 
tions. Actually the system reduces to that of fluid-dynamic stability 
for compressible flows (ref. 15) if we put to = 0 in definition (54) for 
operator L together with the envisaged substitution for q's. This proves 
in the fluid-dynamic space the formal similarity existing between kinetic 
equations (29) and (33). 


IV. GRID-PRODUCED TURBULENCE IN THE LIGHT OF SEPARATION OF VARIABLES 


A. Steady-Flow Formalism Without Isotropy or Taylor's Hypotheses 

Nowadays, a widely used means of understanding experimental data on 
wind tunnel- or grid-produced turbulence is to utilize the theory of 
isotropic turbulence with proper interpretation for obtained solution. 

This method of approach, however, is subject to two assumptions in 
addition to the one for closure. The first of these is that grid-produced 
turbulence is (spherically) isotropic. Speaking more exactly, this is 
an approximation, in nature, rather than an assumption; in fact, a grid- 
produced turbulence is not spherically but cylindrically (or axially) 
isotropic (ref. 16) because the presence of a constant field vector 
u = (U,0,0) (U; the mean velocity in the wind tunnel) distorts the tensor 
field to be so. A consequence of the isotropy assumption is to suppress 
the pressure-velocity correlation identically (ref. 17); 


In the lack of this quantity (and under homogeneity condition) the triple 
correlation is the only term causing effects other than simple trans- 
lation or molecular diffusion of turbulence. (See Eq. (138) of Section 
VI.). The danger is; if a slight anisotropy (in the spherical sense) 
exists in the actual grid flow, causing a slight response in the form of 
nonvanishing r! 1,2) , it may well dominate effects due to the terms of 
0(A 3 ) because it is of 0{A 2 >; a lower order correlation. This means that 
the isotropy postulate needs to be checked carefully before the quality 
of a closure condition to be effected on A 3 is discussed on the basis of 
comparison with experiments done for a grid— produced turbulence. 

The second assumption is Taylor's assumption claiming that the turbu 
lence fluctuation is carried with the mean flow; 
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t 


X 

U 


(57) 


where U is the (constant) velocity of the wind-tunnel flow. This assump- 
tion is needed to interpret space (x ) -dependent experimental data in 
terms of time-dependent solution of the isotropy theory. This assumption 
has no precarious effects as the isotropy assumption does. It would, 
however, be avoidable to invoke if a more direct formalism is available 
such that decay of turbulence is described in terms of space coordinate x 
instead of time t. We will show in what follows how these hypotheses are 
eliminated, and will see what are hidden behind. 


B. A Hidden Solution 

With a frame of reference fixed to a grid or turbulence generating 
apparatus, Eqs . (51) and (52) under uniform (u = constant) and steady 

(3/3t = 0) flow conditions read 


~ (1) 

La' 01 + ^ 




9X 1 
_.j 


= o 


+ I 9 


3q 


(1) 


Bx 


3 9Xj 3x r 


= 0 


(58) 


(59) 


with 


L = 


ioj + U 


3x, 


(60) 


where v denotes the kinematic viscosity. Taking divergence of Eq . (59) 

( 2 ) 

and eliminating the pressure p 


L - V 


3x 


from the two equations, we have 
( 1 ) 


,2 „ (i 

a q 




ax; 


3x . 3x 
3 r 


- 0 


(61) 


J 


or for incompressible flows (q ^ = 0) 


9q (1) 

23 

3x . 

3 


= 0 


(62) 


. iuj + U _ . v 


3x 4 


aV 15 

3 

- 2 

3x, 


= 0 


(63) 
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Note that Eq. (63) is the Orr-Sommer f eld equation under a uniform flow 
condition. The solution of Eq . (63) is written in the form: 



d 0, j 


j 


(64) 


where q Q ^ an ^ q L/ j obey, respectively, a generalized Oseen equation and 
Laplace's equation; 


- iw + U - v — 


9x-. 


3x 


0,3 


= 0 


8 2 q 


L, j = 


ax' 


(65) 


( 66 ) 


Before proceeding further it would be relevant to see at this point 
how Eq. (63) compares with the classical counterpart under the same 
approximation = 0: In the classical theory this case has been called 

the final stage of decay, and the governing equation is 


a 

at 


- 2 v 


R 


Sr: 


( 1 , 1 ) 

j ,1 


(67) 


where r^ denotes x^ - x^. Equation (67) is apparently a diffusion 
equation, parabolic and of the second order, whereas, in contrast, our 
equation (63) is elliptic and of the fourth order. This difference is 
mainly ascribed to the a priori lack of the pressure-velocity correlation 
(56) that is an immediate consequence of the isotropy postulate, causing 
thereby a difference by the factor of Laplacean. 


C. Dispersion Relationships 


The two types of mode given by q . and q . of (65) and (66) , 

O f 3 L ,3 

respectively, are attached simple physical interpretations through dis- 
persion relationships; namely, relationships between the frequency w and 
the wave number k^ . To obtain this we assume q in the form 


= exp ik.x. 

3 3 


( 68 ) 


where k 2 and k^ are taken to be real reflecting invariance of correlations 
with respect to Galileian transformation in ( X 2 /X 3^ plane normal to the 
flow direction. Then for the Oseen mode q^, we have an algebraic version 
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0 


(69) 


of differential equation (65) in view of (68) : 

— i to + i(k 1 ) Q u + v[(k x )Q + k 2 ] = 


with 


2 .2 ,2 

k = ^2 + ^3 


Expressions (69) determines complex wave number 

(k, ) ^ = (k, ) _ _ + i (k 


^l y 0 1 0 , r 


1 0, i 


in the form; 


(k l } 0,r 


= u 


< k l> 0 ,i - 5 “Vo.r + ^ 


(70) 


(71) 


(72) 

(73) 


Since the quantity on the left-hand side of (72) means the wave-propaga- 
tion velocity, this relation shows that the wave propagates at the fluid 
velocity; namely, the turbulence is frozen to the fluid. In other words, 
this mode obeys Taylor's hypothesis. Although (72) as the solution of 

Eq. (69) is subject to an approximation, correction terms to appear on 

2 

the right-hand side is shown to be of 0(v ), proving that this mode obeys 
the Taylor rule exactly within fluid-dynamic regime. Condition (73) , on 
the other hand, describes how this mode decays spatially. 


In a similar way we are led to the dispersion relationship for the 
Laplace mode: Substituting (68) into (66) we obtain 

(74) 

(75) 


(k l>L,r " ° 


(k l } L,i = K 

claiming that this mode is stationary in space, namely, relative to the 
grid, and spatial rate of decay is k. It shows us that this mode decays 
faster than the previous mode (73) (~0(v)), so that it is observed 

immediately downstream of the grid, then is followed by the Oseen mode, 
corresponding to the final stage of decay. Thus we have learned that the 
turbulent field consists of superposition of two modes, the one obeying 
and the other disobeying the Taylor hypothesis, respectively, and that 
the hypothesis is met better as we go downstream because of faster decay 
of the Laplace mode. 
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D. Determination of Axially-Isotropic Correlation Tensor with Given 
Data at an Initial Plane through a Bessel Inversion 

As contrast to the classical theory which is time-dependent and 
spherically isotropic, we aim at formulating the velocity correlation 
tensor as space-dependent and axially-isotropic . For this purpose it 
seems more convenient to write our governing equations (62) and (63) in 
the form 


with 


■ 9 

1 

n 2] 

n 2 (1) 


f x “ + 93Z " 

R 

V 

V q l 

= 0 

• 4- 9 

1 

v 2 
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9q 
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8q 2 
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9z 


= 0 


0 


(76) 


(77) 


(78) 


(x, y , z ) — ( x ^ ' x 2 / x^) 

V 2 = 9 2 /9x 2 + 9 2 /3y 2 + 9 2 /3 z 2 

l 

R = U - 

V 

-> w, (redefinition) 


'N 


) (79) 


J 


where R is the Reynolds number, and the grid mesh size l is used as a 
characteristic length for nondimensionalization. If the solution is 
sought in form (68) , where we have redefined 


Ik. -> k. 

1 1 

Eqs. (77) and (78) are solved for q^^ an< ^ ^3^ as 


( 1 ) ■ -2 

q 2 = 


9q 


(1) 


k 3 q + k 2 ~9F 


q^ 1} = Ik" 


9q 


(1) 


" k 2 q + k 3 “93T 




) (80) 




where q is an auxiliary function obeying eq. (65). We can, therefore, 
put 
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q = B 0 q 0 (81) 

with complex constant B . As for qj 1 ^ we have already obtained a general 
solution 


q 


(i) 

1 


A 0 q 0 + A L q L 


( 82 ) 


where q and q obey (65) and (66), respectively. 

O ti 

Our task here is to show the following: the complete knowledge of 

the correlation tensor 


R U) (X,X) = q (1) (x)qj 1} * {i)dwdk 2 dk- 

3 / -t J 3 z ~ 


(83) 


at any set of points (x,x) is provided if a complete knowledge of the 
tensor element is given at the initial plane 


x = x = 0 

as functions of the transverse separation 


r = [ (y - y) 2 + (z-z) 2 ] 


1/2 


(84) 


and in the form of frequency-analyzed correlation 


( 1 , 1 ) 

j 


<r< 1 ; 1) > d, 

j r t cj 


(85) 


as appearing on the right-hand side of the integral. In mathematical 
language, six independent initial data of the symmetric tensor 
[ R f ^ ^ ^ ( r ) ] in the transverse plane x = x = 0 are necessary and suffi- 

J r t OJ 

cient for unique determination of six constants; A , A , and B Q in (81) 
and (82) . 


The way in which these complex constants are calculated in terms of 
the initial data is described in some detail in Appendix A. A brief 
sketch of the solution through a functional inversion and the meaning 
of the obtained results are discussed here: From Eqs. (83), (68), and 

(82) we have 




2tt 


ik-. x 

(A 0 e 


n -<x w * - ik l* = 
a l g > (A o e 


+ A^*e KX ) Jq ( xr) KdwdK 


( 86 ) 
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where k^ is given by (71), denotes Bessel's function, and an integra- 
tion with respect to the azimuth angle arctan k^/k^ has been carried out. 
At the initial plane x = x = 0, and in terms of the frequency-analyzed 
initial data of (85) , Eq. (86) reads 


[R 


( 1 , 1 ) 

1,1 


2tt 


% + a l> (A J + A £ )J o <Kr,KdK 

0 


This relationship is inverted by means of an integral formula 
give 

CO 


(A o + V (A o + a l } 


2tt 


[R 


( 1 , 1 ) 

1,1 


w rJ 0 (<r)dr 


(A. 9) to 


(87) 


Six independent relationships of similar forms are obtained likewise. 

These determine constants A , A , and B uniquely depending on two para- 

Ob O 

meters; the frequency u) and the transverse wave number k. 

Approach along this line, in principle, to the grid-produced turbu- 
lent phenomena has been initiated by Nakagawa (ref. 18). The formulation 
described here is essentially the same except that in the earlier work 

the Laplace mode was suppressed (A - 0) , so that the initial-value 

b 

problem with the full tensor components as posed here would have caused 
overdeterminacy . 


E. Pressure-Velocity Correlation as a Prime Order Quantity 

Once the amplitudes of the two modes have been determined uniquely 
((A. 10) through (A. 15)), we can show that the pressure-velocity correla- 
tion 


r! 1 - 2 ’ 

J 


q_f^(x)p^ (x)doo dk^dk^ 


( 88 ) 


is necessarily nonvanishing and behaves as a correlation of prime order 

(12) 

of magnitude as contrast to R; ' = 0 of the isotropic turbulence, alsc 

(12) ^ 

to R; ' - O(A^) of some other cases of homogeneous turbulence. 


( 1 ) 


Solutions (82) and (80) for q. 

( 2 ? 

separated pressure fluctuation p give 


and Eq. (59) solved for the 


[ 2 ) 


= -(1 + iu)/K)A q 


(89) 
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An important implication of this formula is that only the Laplace mode 

contributes to the pressure fluctuation, namely to generation of noise, 

so that the Oseen mode is interpreted as representing quiescent component 

f 2 ) 

independent of its turbulence intensity. Since the quantity p is 
proportional to noise amplitude human timpanic membranes perceive, this 
observation serves to give an intuitive criterion as to whether the 
Laplace mode is important under an experimental condition: It would be 

non-negligible compared with the Oseen mode if noise is sensible in the 
process of generation of a grid-turbulence. According to this criterion, 
the Laplace mode seems to play a role even for relatively low-speed flows 
through a grid such as an outlet flow of a room aircondi tioner . 


In view of expression (89), the pressure-velocity correlation (88) 

2 

is seen to consist of terms of 0 (A A ) , 0 (A ) , and 0(B^A T ). Since none 

U L L 0 L 

of these terms vanish other than exceptional cases of no particular 

interest we may conclude that the pressure-velocity correlation has the 

same order of magnitude as other double correlations. This conclusion, 

f 1 2 ) 

as it contrasts drastically with the formula Rj ' } ~ O(A^) used for 

homogeneous turbulence, will need a supplementary explanation to reconcile 

the seeming gap: If a flow is homogeneous but not spherically isotropic, 

( 12 ) 

R. ' is shown to obey Poisson's equation of the form 


9x 


i -J 1 ' 2 ' - 0<* 3 > 


where the source term comprises only the triple correlations. Its general 
solution is given, accordingly, by 

r ( 1/2) _ p + H 
D 

where P(~0(A^)) is a particular solution of the Poisson equation, and H 

is a harmonic function to be responsible for the boundary conditions. It 

(1 2 ) 

should be borne in mind that the assertion R^ ' - O(A^) above holds 

only under the condition H ~ 0 or O(A^), which is the case with the 
homogeneous turbulence with boundary conditions given at infinity. The 
fact that we are dealing with a flow which is semi-infinite (x J> 0) in 
extent rules out the possibility, because the function H should cope with 
a boundary condition that is imposed at a finite point (x = 0) , and is 
of 0 (^ 2 ^ as has been assured. 
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V. APPLICATION TO INCOMPRESSIBLE WALL TURBULENCE 
A. The Reynolds Stress as Coupling Function 

In this section we will focus on incompressible shear turbulence 
with special reference to wall turbulence, and will see how the equations 
governing fluctuations are coupled with those describing the mean motions 
by demonstrating the actual expression of the coupling, namely, the 
Reynolds stress. 

Under the condition of incompressibility (p: const., = 0) 

equations governing average and fluctuating quantities constitute a 
closed set with Eqs . (42) , (43) , (51) , and (52) ? 


3u /3x = 0 

r r 


- -"if 

3q r /3x r - 0 


( _ iw + gt + u r 9x 


q . q* dw} 
^ j ^r 


9 a 2 x 3u -j ,i 3p <2) 

- v — q. + 7777^- q_ + - — = o 


8* 2 ^ 9x r 
r 


p 3x . 

3 


(90) 

(91) 

(92) 

(93) 


where (NS) j is the Navier-Stokes ' operator (2) replaced with the average 
quantities, and q^ is defined by 

-1 ( 1 ) 

q . = p q 
3 J 

If, in addition, the flow is considered to be a parallel flow 


u . = U6 . , u (x~) 
3 3i z 


(94) 


with X 2 **axis taken normal to the wall surface, we can assume, without 
loss of generality, as 


= Yj(y)exp i(ax + Yz) 


where 


(x , y , z ) = (x 1 ,x 2 ,x 3 )/£ 


(95) 


(96) 


denotes a cartesian coordinate nondimensionalized by a fluid-dynamic 
characteristic length £ . In expression (95), a and Y are wave numbers 
in the directions of the velocity (x) and the (minus) vorticity (z) , 
respectively , made nondimensional in terms of £. 
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Since, in our system, the quantities u ^ and are dependent on each 
other through Eq . (91), the parallel flow assumption (94) prohibits the 

Reynolds stress to appear on its right-hand side from depending on x 
and/or z explicitly. It leads necessarily to the requirements that both 
a and Y be real. In fact, if a is a complex number, the Reynolds stress 
would be x-dependent like exp [-2 1 (a) x ] , where I( ) stands for the 
imaginary part. Of the two requirements, reality of Y is understood as 
a postulate for invariance of the Reynolds stress to a Galileian trans- 
formation in the y-z plane. Reality of u, on the other hand, limits 
permissible classes of solutions to those u(y)'s which have neutral 
stability under a certain condition. It will be discussed in Section V.D. 


Provided, for the time being, that the condition is fulfilled, Eqs. 
(90) through (93) reduce to a system of ordinary differential equations 
in y in the following fashion: First we introduce two auxiliary variables 

Y and Z to represent Y_. of (95) obeying constraint (92) of the continuity 
condition ; 


Y 


1 


_i 

a 


(DY + 



Z) 


Y = Y 
1 2 1 


Y 3 " " 


ie 


a (1+0 ) 


(97) 


where we have introduced the notations 

D = d/dy 

and 

0 = y/a (98) 

the latter representing a three dimensional effect of the fluctuations. 
Equations governing Y and Z are, then, obtained from (93) by eliminating 
the pressure fluctuations, yielding (ref. 10) 

[ (u-c) (D 2 -A 2 ) - jig (D 2 -A 2 ) 2 - u " ] Y = 0 (99) 

y 

[u-c - (D 2 -A 2 )]Z = -A 2 I [u-c - -r-+r (D 2 -A 2 )]Y dy (100) 

XOtK J lCtR 

where u' stands for Du, and 
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( 101 ) 


R 

A 


U l/v 


(1 + 


eV /2 


a 


c = w/a 


respectively, give the Reynolds number, the magnitude of the wave number 
vector, the streamwise propagation velocity of the coherent wave as 
discussed in Section IV. C. In the latter expression w is the nondimen- 
sional frequency subject to redefinition (79). 


Boundary conditions for Y and Z are obtained through the fact that 
the fluctuation amplitude (97) must vanish at the solid boundary or at 
an asymptotic limit where the flow tends to equilibrium; 


Y = DY = 0 
Z = 0 


( 102 ) 


Equations (99) and (100) for Y and Z constitute a closed set together 
with the Navier-Stokes equation (91), namely, 

(NS) X = -dR 12 /dy (103) 

where 


r 12 (y) = 


-1 G * 

a (DY + 9 Z) Y dw} 

1+0 


(104) 


is the shearing component of the Reynolds stress nondimensionalized 

properly. Without solving for Y and Z, it is assured from condition 
2 

(102) that Y~y and Z~y for y<<l, so that the Reynolds stress (104) 
varies like y^ in agreement with the known dependence (ref. 19) in the 
viscous sublayer. 


B. A New Scaling of the Orr-Sommerf eld Equation 
Appropriate for Wall Turbulence 

The fact that the amplitude Y for the y-component of the turbulent 
fluctuation obeys the Orr-Sommerf eld equation as evidenced by Eq . (99) 

has been pointed out in reference 2. For a fully-developed wall turbu- 
lence, however, the equation has a peculiar difficulty in the method 
of solution caused by the fact that the curvature term u" of the equation 
has an extraordinary magnitude, not to be neglected as is disposed of 
through the classical scaling (refs. 20,21). A key to the solution is 
provided (ref. 10) by employing a new scaling parameter e in stretching 
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coordinate from the physical (y) to the scaled (n) ones; 


n = y/e 


e = 


(aR) 


- 1/2 


( 105 ) 


-1/3 

Note that the classical scaling is e - R , so that the rate of 
stretching is greater for the present case. With this scaling law the 
Orr-Sommerfeld equation (99) reduces to 


Y- i ( u-c ) Y + iuY = 0 (106) 

with 

Y = dy/dn f etc. 

retaining the curvature term as a leading one, and in a form as if added 
to the two-term equation of Heisenberg (ref. 20) widely used in the 
stability theory. 

Because of the non-negligible third term intervening in the equation 
the classical asymptotic method which solves the two-term equation using 
universal (Hankel) functions is no longer available here. It turns out, 
however, that an alternative asymptotic method is workable through an 
observation that owing to the presence of the third term Eq. (106) is 
integrable once, yielding 


Y* - i (u-c) Y + iuY = iC ]L (107) 

Four independent solutions <l>^ of this equation are obtained with some 
similarity to and more essential difference from the corresponding 
solutions of the classical equation. For the method of solution and 
their actual expressions see Appendix B. A few remarks are in order: 
Solutions and $ 2 ' the counterparts of slowly varying solutions of the 
classical equation, behave actually so in and only in asymptotic limit 
(n >> 1) , 


~ - (u-c) 


(u-c) 2 dn 


(108) 


$ 2 ~ (u-c) 


In the region of high shearing rate (n ~ 0(1)), these are subject to 
rapid variation in contrast to 'uniformly' slow variation of the 
classical solutions throughout the field. Two other solutions and $ 
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represent, respectively, rapidly decaying and growing functions with the 
asymptotic expressions in agreement with those of classical solutions 
given in terms of Hankel functions of first and second kinds. The 
solution procedures are: functions $ ^ and $2 are obtained by solving 

first order equations of which the critical layer (u-c =f 0) is a regular 
point. No multiplevaluedness as encountered in the classical theory 
occurs in the quadrature through the point. Once these functions are 
determined through a straightforward numerical method, and are 
solved analytically and exactly in terms of them. 


C. An Eigenvalue Problem 


A salient feature of our scaling law (105) is that the scaled or 
'inner* equation (106) serves also as the starting equation in a system- 
atic method of successive approximation solving the full Orr-Sommerf eld 
equation (99) . In fact, expansion of the eigenfunction in the following 
form; 


y = y 


( 0 ) 


+ (eA) 2 Y (1) 


(eA) 4 Y (2) 


y = en 


(109) 


e 


(aR) 


- 1/2 


substituted in the Orr-Sommerf eld equation yields the following set of 
equations determining successively? 


L (0) Y (0) - 0 


L (0) Y (1) = L (l) y (0) 


L (0) Y (n) = L (1) y (n-1) _ Y (n-2) 


(n i 2) 


( 110 ) 

( 111 ) 

( 112 ) 


where operators and are defined by 

= d 4 /d 4 - i(u-c)d 2 /dn 2 + iii (113) 

L (1) = 2d 2 /dn 2 - i (u-c) (114) 


Equation (110) , together with (113) , confirms the envisaged assertion not 
to be expected for its classical equivalent. This property enables us 
to reach the full solution F simply through advancing the successive 
approximation starting with the inner solution 0^ : 

F = $ + y ( eA) 2n 4> (n) ( i :1 , 2 , 3 , 4 ) (115) 

1 1 n=l 1 
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The higher order approximation is obtained by a method of finding a 

particular solution of Eq . (ill) or (112) as given in Appendix B. 
Asymptotic expressions (n>>l) of and F 2 thus obtained are of particular 
interest f 


F 


1 


- (u-c) 


Y 

(u-c) 2 dy + A 2 


y 

(u-c) 2 dy * 

' A 

y 


y y 

(u-c) 2 dy (u-c) 2 dy + A 4 . .. 

/V ^ A 

y y 


F 


2 


(u-c) 


1 + A 


2 


y 


dy 

(u-c) 2 


y 

2 4 

(u-c) dy + A 

y 


(for y>>e) 


These expressions which follow immediately from Eqs. (108) and (B.37) of 
Appendix prove exact coincidence with Heisenberg's outer (inviscid) 
solution (ref. 20). This fact indicates that F of (115) thus obtained 
is valid for the outer- as well as the inner-layers, enabling us to by- 
pass the conventional procedure of matching the two solutions. 

With these F's having so derived the eigenfunction Y is written as 

4 

y = y C F 
x = l 1 1 

Since, as seen in Section V.B, F^ grows drastically with r\ , we need to 
put = 0 to cope with boundary condition (102) at n = °°. On the other 

hand, condition Y(0) =0 imposes = 0, because F^tO) = 1^0, whereas 

F 1 (0) = F 2 (0) =0. Thus we are led to the final expression of the eigen- 
function ; 

Y = C^F + C 2 F 2 (116) 

depending only on two functions and F 2 throughout the region in 
question. Note, for comparison, that its equivalent of the classical 
asymptotic theory comprises three terms (C^ ^ 0) . Imposition of two 
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remaining conditions 


Y' (0) = 0 


(117) 


Y' + AY = 0, (at y = y) 


(118) 


yet to be satisfied by Y of (116) leads to an eigenvalue condition in 
the following form, 


(0) 

F'(0) 

(y) + AF 1 (y) 

Fj(y) + AF 2 (y) 


(119) 


Adequacy of using the outer boundary condition of form (118) in turbulent 
cases is discussed in Appendix B. 
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VI. THE SECOND METHOD OF SEPARATION OF VARIABLES: 

NONLINEAR CONVOLUTION 

In this section we will seek a condition which leaves the triple 
correlation, as it appears in the double correlation equation, untouched 
and under which the double correlation equation is still accessible to 
the method of separating variables. 

It has been shown in reference 28 that Eqs. (13) and (14) which, 
according to the BBGKY formalism, are derived from Liouville's equation 
through multiple integrations are also derivable on the basis of the 
following equation; 

B((0 = -£r + V -- 6 (z) -J(z|z)U(zH(z)] = 0 (120) 

3 xJ 

This equation is the Boltzmann equation in which the Boltzmann function f 
is replaced by the microscopic density namely, its fluctuating equiv- 
alent. This equation may well be called the ’master 1 Boltzmann equation 
in the sense that the whole system of the BBGKY hierarchy is derived also 
from this equation through a more straightforward procedure. In fact, it 
is easily confirmed that Eq. (13) is expressed in terms of Eq . (120) as 

<B U ) > = 0 (121) 

a 

with = /$(z^). Similarly, the two-particle equation (14) with the 
triple correlation retrieved is written as 


<A 6b B(i a ) + A ^a B(i b }> = <A <Sb“a (A|5 a ) + ^aV^b^ = 0 


( 122 ) 


with operator (2 defined by 


V A «a> E [jt * ? a ' J < 2 a l : 2 a> «a 4 *a + * A <a S <a> < 123 > 


Obviously ft differs from ft of (18) by the presence of terms quadratic in 
a 

A((. Because of this term Eq. (122) includes triple correlation of the 

form Af! k(,~bL, therefore recourse is made to the equation of the 
a a b 

succeeding hierarchy for its determination. This equation is given by 

< l A ( 5 b A^ c B( 1 $ a )> = 0 

a , b , c : 

cycl . permut. 
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( 124 ) 


or, alternatively, in terms of operator Q of (123); 

cl 


l 

a , b, c : 

cycl.permut. 



<Ai5 b A)5 c >J(z a |z s ) (<A^ a A(S 



0 


where the summation is over cyclic permutation (a+ b-> c-> a) , each term 
yielding three of them. Equation (124) comprises three kinds of corela- 
tion quantities; triple and quadruple correlations in linear forms and 
products of double correlations. 


Let us search for the solution of Eq. (122) again in the form of 
separation of variables 


< H a H h > 


1 (z ;ta ) 4> (z, ; w, ) 6 (w + a>,)dw duv 
aa b b a bab 


(125) 


which is an alternative expression of (16) with supplementary condition 
(20) . In the same way, we assume the triple correlation governed by 
(124) to be separable. 


<A ^ a ^b^c > = 


6 (Eto ) n [<t> ' (z ;<d ) dm ] 
a a a a 


(126) 


where the sum (Z) and the product (II) are over three indices a, b, and c. 
An integral expression of the like is assumed for the quadruple correla- 
tion with the integrand given by a quadruple product of + In general, 

4- a 

a function 0 of (125) will differ from 0 of (126 ) r etc. 

Consider an iteration procedure in which the product of double 
correlation, namely, the last term of (124) is neglected at the outset, 
whereas no approximation is invoked in (122) . Then, it is readily seen 
that Eqs. (122) and (124) reduce to the same equation if we put 

<t> = 4> + = <f’ + + (127) 


with the resulting equation 


-irn + ^ + v 


_a_ 

-y 

dx 


- j(zjz)[f0 4- f 0 + 0 (aj- 0 ) 1 ) 0 (u) f )du>* ] = 0 (128) 


where we have dropped the subscripts (a,b,c) because there would occur 
no confusion from this point on by so doing. If the nonlinear term 
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(with the integral) is deleted in Eq . (128) in line with the ternary 

chaos hypothesis, it reduces to Eq . (IB) as it should. 

Looking back the procedure leading to the closed Eq . (128) for <j> , 

we see that the second equality of condition (127) serves as the closure 
condition. In fact, only under this condition, the governing equations 
for $ and prove to be identical, both yielding Eq. (128) , therefore 
the first equality of (127) is justified a posteriori. Since it is 
intuitively reasonable to interpret mathematical realization of a physical 
concept of taking average < > by the integral of form (125), (126), etc. 

then our closure assumption that triple and quadruple correlations are 

t 1 1 

expressed using a common variable <J> = <)> seems to have a sound basis. 

Physical implications of the nonlinear term may be easier to under- 
stand by looking at the fluid-dynamic version of Eq . (128). The moment 

expansion to lead to fluid equations was described in Section III, so 
will not be repeated. Only the final forms of the equations governing 
the correlation tensor for incompressible flow are given here: 


(x,x) = 

j£ 


q . (x; u) g (x; -u) dco 
1 t 


: 129 ) 


3q 


3x 


2 = 


-iu) + ttt + U 


3t r 3x^ 


- v 


3x 


3u . 

‘Jj + a^T q r 

J r* 


(130) 


1 3jd 


( 2 ) 


P 3x 


3x 


q . (u)-u)’ )q (u)’ )du)' 
’ j r 


= 0 


(131) 


where a. has been defined in Section V.A, the q. 's except those inside 
3 3 

the integral of (131) stand for g. (x;w). This equation is to be con- 
trasted with Eq. (93) where the triple correlation has been ignored. 

For a closer look at Eq. (131) we assume the mean flow to be a 

2 2 

parallel shear flow [u. = S^utx^], and apply an operator 6 2 9 - 

3 2 /3Xj3x 2 to the equation to eliminate the pressure term. We have, then. 


38 



r 3q 2 (ww* ) 


q j (oj ' ) dw 


Aq 2 



3x . 

3 

3q r ((jj-U)' ) 

3x . 

J 



doj 1 


( 132 ) 


where the operator A on the left-hand side is given by 


A 


ioj + 






( 133 ) 


that is, the Orr-Sommerf eld operator with added term of -iw reflecting 
the effects of the variable separation. Thus we see that the velocity 
variable q 2 obeys a (generalized) Orr-Sommerf eld equation, 'driven' by 
terms due to the triple correlation that are quadratic in fluctuation 
velocities with different internal variable, namely, different frequencies. 
Thus, the role of the triple correlation has turned out to be the following: 
For not too weak, turbulence in which the terms of 0(q) still stand, it acts 
as controlling the interaction between turbulent fluctuations with dif- 
ferent frequencies; the form of the nonlinear interaction is mathematically 
a convolution integral and physically a continuous cascade of turbulent 
energy in the frequency space. This is in contrast with the role of non- 
linear terms of the Orr-Sommerf eld equation of classical stability theory, 
where the energy cascading is through discrete higher harmonics (refs. 11,12). 

VII. CONCLUDING REMARKS 

In concluding this report, a brief review from a different angle 
will be given to help clarify physical implications of the two closure 
assumptions we have employed; the neglected triple correlation (Section II) 
and the identity of separated variables in triple and quadruple correla- 
tions (Section VI) . One way of doing this would be through Eq . (1) for 

incompressible flows written in the following form: 
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(4) 


We note that this equation is exact within this context, but is not 
closed because of the triple correlations 

0 (1,1,1) ,+ ~ i . _ 3 , 

Rj r ^ (x,x,x) = p <AUj Au^Au^> 


intervening in the form 


( 1 , 1 , 1 ) 

j r,l 


( 1 , 1 , 1 ) 

j rl 



also R. n defined likewise, depending on two space points. Equation 

3 ' r<c (11) 

(134) describes how the double correlation R^' / in the course of con- 

vection and diffusion (terms (0) and (1) , respectively) , interacts with 
the mean shear flow (2), with the pressure-velocity correlation (3), and 
with the triple correlation (4) . 


For the classical isotropic turbulence terms (2) and (3) vanish; 
thus, the triple correlation is the only factor responsible for active 
mechanism in turbulence evolution as contrast to passive ones through 
the convection and the diffusion prevailing in the final stage of decay. 
Grid-produced turbulence as formulated in Section IV., in contrast, keeps 
another active term (3) finite, lessening the danger of making assumption 
on term (4), since the latter term is higher order in correlations. Its 
finiteness leads to finding of a pre-final stage where turbulence decays 
more quickly than downstream (the Laplace mode) , within the same regime 
of ignored triple correlation. 
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For shear turbulent flows where terms (2) are also finite, the 
danger of neglecting the triple correlation is diluted more than in the 
grid- turbulence case. 

The attempt to incorporate effects due to finiteness of the triple 
correlation through the second assumption is of exploratory nature. 
Term-by-term comparison of Eq . (134) with Eq . (131) tells us that term 

(4) is represented by an integral with respect to the frequency in the 
form of a nonlinear convolution, describing a peculiar mechanism of 
direct cascading of energy from one turbulent component to the other with 
the frequency arbitrarily far apart. This is to be contrasted with in- 
direct ones due to term (2) , where the interaction between two modes 
are via a mediator, namely, the shearing motion of the mean flow. 
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APPENDIX A 


UNIQUE DETERMINATION OF COMPLEX AMPLITUDES A 0 , A L , AND 3 0 OF 
GRID-TURBULENCE FLUCTUATIONS IN TERMS OF CORRELATION DATA 

AT AN INITIAL PLANE 


The problem to be solved separately here is to determine complex 

amplitudes A , A , and B appearing in solutions (81) and (82) of the 
O L O Q\ 

separated fluctuation variables q . in terms of a set of correlation 
data given at an initial plane; 


x = x = 0 


(A. 1) 


( 1 ) 


To work out this task we first write q^ given by (82) and (80) in 
the form 


( 1 ) 




q i = {A 0 X 0 + A L X L )SXP iU 2 y + k 3 Z) 


( 1 ) 


-t [ (k 3 B Q + ik 2 k 1 A Q )X 0 - Kk 2 A L X L ]exp i (k 2 y + k 3 z) \ (A. 2) 


<i) _ _L 

q 3 2 

K 


[ (~ k 2 B o + ik 3 k l A o^ X 0 ~ <k 3 A L X L ]exp i (k 2 y + k 3 z) 


with 


X (x) = expUk-^x) 

X (x) = exp (- kx) 

L 


(A. 3) 


where we have employed Eqs . (81) , (68) , (72) through (75) . With these 

formulae expression (83) with (85) for frequency-analyzed velocity 
correlation tensor is 


[R. J = q, ; (X) q; ±; (x)KdKdB 

J t f u) J J e 


( A . 4 ) 


where we have suppressed superscript (1,1) in Eq. (83) and 6 is defined 
by 


k 2 ~ k cos B 
k^ = k sin 3 


(A. 5) 


We introduce a cylindrical coordinate system in the fluid-dynamic space; 

y - y = r cos 0 
z - z = r sin 0 


(A. 6) 
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Since we may assume constants A ^ , A^, and B Q not to depend on 0 because 
of cylindrical isotropy of the phenomena, integration with respect to 
variable 0 - 0 in (A. 4) is seen to be practicable: To carry out the 

integration we make use of the following formulae; 


2tt 

► 

0 



[iicr cos 


(0 - 0)]d(0-9) 


2TriK 


cos 

sin 



(<r) 


4 



[iicr cos (0-f3) ]d ($-0) 

= 7T < 2 [J 0 (<r) + cos 20 J 2 (<r) ] 


2tt 

f 2 

k 2 k^ exp [ ixr cos (0-8) ]d (8-0) = -irk sin 20 J 2 (<r) 

0 


'N 


> (A. 7) 


J 


where J's are the Bessel functions and r has been defined by (84). In 
obtaining these formulae we have made use of Hansen's integral representa- 
tion for Bessel's function. 


J n (x) 



e 


ix cos0 


cos n0 d9 


Assume that we have a set of measured data for six (independent) 
tensor components at x = x = 0 and 0=0. These measured data 

given as functions of r are shown to be necessary and sufficient for 
unique determination of the complex amplitudes in the following way: 
Substituting (A. 2) into (A. 4) with integral formulae (A. 7) into account, 
we have the following expressions: 
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oo 


R 


llw 


= 2tt 


(A + a l ) (A* + A*) J 0 (Kr}KdK, 


R 


22w 


- 7T 


-1 


k d k { B q B* [ J q (< r) + J 2 (Kr)l + k 1 k 1 *A Q A* [ J Q ( Kr ) 


- J 2 (<r)] + iKk* A ! A* [ J Q (kit) - J 2 (<r)] - iKk^A Q A* [ ( <r ) 

- J 2 (<r) ] + k2a l a l [Jq (< r) - J 2 (<r) ] }, 


R 


12u) 


= ~2tt 


R 


13oo 


— 2 7T | 


dK J 1 (Kr) (A Q + A l ) ( ik*A* + KA*) , 


dK J ]L (Kr) (A Q + A l ) B* , 


) ( A . 8 ) 


R 


3 3oj 


= Tl 


k 1 dK{B Q B* [ J Q (<r) - J 2 (<r)] + k i k i A o A o ^ J 0 + J 2 ^ <r ^ 


+ iK(k*A A* - k A*A ) [J (icr) + JjKr) ] + k A T A* [ J (<r) 

1LO 1 L O U 2 LLo 


+ J 2 (<r) ] } , 


R 


23(d 


k *dic{B (ik*A* + kA*) + B^-ik^A^ + kA^) } J 2 (iar) 


J 


On the right-hand side of these expressions, the real part should be 
taken because of the reason mentioned in Section II. If we are reminded 

of the fact that R . 0 is the only functions of r, also that the variable 

j to) 

r appears only through the Bessel functions on the right-hand side, we 
see that these equations can be inverted through the formula 


oo 


0 


J _ (ar) J (br)dr 
n-1 n 


r .- 1 




(for b > a > 0) 

(for a = b) 

(for a > b > 0) 


> (A. 9) 


J 


44 



In fact, multiplying the first equation of (A. 8) with J^br) integrating 
with respect to r, and applying (A. 9) we have 


(A o + V (A o + a l ) 




R ll» rJ 0 (Kr,ar 


(A. 10) 


In a similar way the following equations are obtained: 


K 1(A 0 + V (ik i A o + KA ^ = ~ h 


R 12 u) rJ l (Kr)dr 


(A. 11) 


-< _1 < A 0 + V B S = h 


R, _ rJ. (icr) dr 
1 / 0 ) 1 


(A. 12] 


- k iVS> + K( VS + vi> = 


2k 

7T 


R 23u) J l (Kr)dr + ~ 


R 23m rJ 0 (Kr)dr 


(A. 13) 


B 0 B 0 4 TT 


0 

2 
K 

4 7 T 


(R 22 W “ R 33J l ~ 2KJ l (Kr) + K rJ 0 (Kr)]dr 


(R 00 +Rn-, ) r J ( kt ) dr 
22 co 33 uj 0 


(A. 14) 


{-ik 1 A 0 + kA l ) ( ik*A^ + kA*) 


4 TT 


(R 22 o)" R 33 o)^ t 2 K J x (<r) 


k rJ^ (Kr) ]dr 


— f 

4 tt j 


+ ifif (R 22 «D +R 33 U )rJ 0 (Kr)dr 


(A. 15) 


where, as before, the real part should be taken on the left hand side. 
We are, thus, led to six relationships (A. 10) through (A. 15) for unique 
determination of the complex quantities A^, A^, and depending on two 
parameters co and k and in terms of experimental values of given as 

a function of r. 
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APPENDIX B 


A NEW ASYMPTOTIC METHOD OF SOLVING ORR-SOMMERFELD EQUATION 


This appendix is mainly concerned with obtaining four independent 


solutions $ ( i : 1 , 2 , 3 , 4 ) for Y; 


y = c 1$1 + c 2 * 2 + c 3 $ 3 + c 4 o 4 


(B.l) 


obeying Eg. (107), namely, 


Y - i(u - c)Y + iuY = iC 1 


(B . 2 ) 


For the procedure of deriving this eguation from the full— term Orr— 
Sommerf eld equation, see the main text (Section V.B) . In numbering the 
four independent solutions, we have followed the conventional usage, for 
example that employed in reference 21, so that each counterpart have the 


and 0 ^ tend to slowly varying functions. 


same asymptotic behavior; uuu TT 2 

whereas $ 3 and <J> 4 exhibit exponential decay and growth, respectively. 

Our procedure of obtaining the solutions will be carried out in the 
order of $> 3 •+■ 4> 2 $ 4 -*■ 4^ in such a form that each solution will be 

worked out using knowledge of 4> ' s at previous stages, and 4> 4 and 4>. will 
be solved analytically in terms of $ 3 and 4> 2 . 

Solution 0^ 

First we will consider homogeneous equation (C^ = 0) of (B.2) which 
provides with solutions ^ through ? 


Y - i { u - c)Y + iuY - 0 


(B. 3) 


This equation has a structure such that the critical layer (u = c) is an 
ordinary point, so that no artifice as needed in the classical treatment 
is necessary in integrating the equation across the point. In view of 
this feature, also of the prospective exponential decay, the solution ^ 
may be assumed in the form 


Then Eq . (B. 3) 

reads 


= exp 


X dn (B.4) 

subject to this transformation of the dependent variable. 
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ir, alternatively and more conveniently , 

A = y - i (u - c) , 

. 3 

y = - 3A y + 4i(u - c) A - A 


(B.5) 


Since we are seeking a solution of (B.5) which decays exponentially in 
the asymptotic limit (A, y ~ 0) , the root for A with negative real part; 


A 



c) 


1/2 


y - i (u - c) 


(for v) >> 1) 


( B . 6 ) 


is seen to meet the purpose. In solving Eq . (B.5) for A and y with 

asymptotic conditions (B.6) it is advisable to introduce a variable- 
transformation 


2 

V = ~y + 2i (u-c) - A 

2 

W = -2y + 3i (u-c) - A 


(B . 7 ) 


and to work with the equations in the new dependent variables (V,W) , 


V + A V — 2iu = 0 

W + 2A W - 3iu - 0 
subject to asymptotic conditions 

V, W ~ 0 , as n + 00 



( B . 8 ) 


(B. 9) 


The actual integration may be started at a point n = n(>> 1) where 
initial values 


V = 2iu/A 
W = 3iu/2A 


(B.10) 


with A given by (B.6) are moderately small. The quadrature marches 
inward from this point on with nonlinearity of Eq. (B.8) appearing in the 
form 


A 



(u-c + 2iV - iW) 


1/2 


(B.ll) 


taken into account for n < n* 
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It should be remarked here that this method of numerical integration 
of original equation (B.5) works successfully only for the decaying 
solution $ 2 ' an< 3 a formally identical procedure for the growing 

solution $4 suffers from drastic numerical instabilities. The difference 
has its origin in the fact that the point n - ^ is a saddle singularity 
of Eq. (B.5) for the case treated, whereas in the case it proves to be 
a nodal singularity from which an infinite number of solutions emerge. 

It can be demonstrated that our solution (B.4) has the same asymp- 
totic form as its classical equivalent for n £ 0 >> 1. In fact, a 
straightforward calculation from (B.4), (B.ll) and (B.10) leads to 


$ 3 (n) 


. \ u-c 

= V n) k=z 


5/4 


exp 


-1 | <u- 


1/2 , 
c) dn 


(B . 12 ) 


n 


(for n > n >> l) 

On the other hand, the classical counterpart of $3 is given by (ref. 20) 


i~4>, (x ) 

= ( dx 

L. 3 J 

Cl. J 


, 1/2 u (1 ) 

dx x H 1/3 


00 00 


2 , . ,3/2 

3 Ua o T) 


(B . 1 3 ) 


where H i ^ denotes Hankel function of the first kind, and t and a 
1/3 o 

defined by 


t = (y - 


a = (u* 
o c 


y c ) (otR) 

1/3 


1/3 


with subscript c signifying the value at the critical layer. For i >> 1, 
(B.13) takes the form 


3 


(t) 


'j 


Cl. 




(a o T) 


3/2 


in qualitative agreement with (B.12). 


(B.14) 
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Solution $2 


Once we have obtained one of the solutions of the third order 
equation (B.2) difficulties are considerably lessened as to solving for 
the rest of them because of the following theorem: If n independent 

solutions of a m-th order linear differential equation are at hand, the 
(n + l)-th solution is obtained through solving an equation of (m-n)th 
order. This theorem applied to the current case for the second solution 
$2 warrants an equation of the order 3 - 1 = 2 to be obtained by putting 




2 


n 

G dn 

4 

0 


(B . 15 ) 


and by substituting into Eq. (B.2). This, in turn, is equivalent to 
claiming that if G is assumed in the form 


n 


G 


exp 


(- A + S)dn 


0 


(B . 16 ) 


the transformed variable S obeys a first order (nonlinear) equation. A 
simple calculation actually confirms the assertion, leaving with a Riccati 
equation 

S+S 2 +XS-W=0 (B . 17 ) 


where A and W have been solved from (B.8) and (B.ll), respectively. Of 
the two asymptotic roots of (B.17) for S ~ 0, the one that vanishes as 
n -*■ 00 / namely 

S = W/A (B.18) 

turns out to be the correct choice. For this root alone gives rise to 
the solution ^ meeting the requirement of slow variation as n + 00 • In 
fact, then 


<f> = A 

^2 3 


dr) exp 


(- A + s)dn 


= b 2 (u-c ) , (for n > n >> D 


( B . 19 ) 
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with 


n 

^2 = 3//2 exp - j i + j S dn 

0 

where use is made of the fact that the integrand in (B.19) is a 
rapidly growing function, so that only the portion in the vicinity of 
the upper bound contributes to the integral, then the integral may be 
replaced with 

n 

-1 f f 

(-A + S) x exp (-A + S) dn 

0 

It seems to be interesting to compare $ 2 with its classical counter- 
parts 

U 9 = 1 + 0 (aR ) _1//3 

1 J Cl.,in 
/ 

$2 = (u-c) + O(A^) 

k * Cl . , out 

representing inner (viscous) and outer (inviscid) solutions due to 
Heisenberg (ref. 20), respectively. We note that the classical $ 2 is 
uniformly slowly varying, whereas < 1> 2 as given here is moderately varying 
for n 0 ( 1 ) (so is rapidly varying as a function of the physical 
coordinate y = en) f and tends only asymptotically (n >> 1 ) to a slowly 
varying function (B.19). Also to be noted is the coincidence of (B.19) 
with the leading term of the classical outer solution, proving the sound 
basis of the asymptotic scheme adopted here. 

Solution $4 

Repeated use of the foregoing theorem assures that the third 
solution satisfies a first order linear differential equation to be 
deduced from (B.2). This assertion can be materialized by means of the 
method of variation of constants; 

<&r=r(!)+r$ (b.2 

4 2 2 3 3 



50 



where $ 2 and ^3 are functions of n to be determined. Having eliminated 
$2 between original equation (A. 2) and the supplementary condition 


^ 2^2 + r 3 lI> 3 


0 


standard of this method, we are led to an equation for that is 
essentially of the first order; 



*2*3 " V2 


V 3 " V 3 


This equation is analytically integrated to give 

n n 


4 y 2 


$ 3 W dn + <S> 3 j 
0 C 

where W is the Wronskian formed by $ 2 and 


$ 2 W ^ dr) 


(B . 22) 


W 



(B . 23) 


It is easily confirmed from analytical solution (A. 22) for that 
it has an asymptotic form of exponential growth with 1 % In fact, using 
the same approximation as has been utilized in deriving (B.19) we have 
for (A. 22), 


with 


$ 


4 ! *S*3 


- (u-c) ^^exp 


TT . 

e 


(u-c) 1/2 dn 


9 






(B . 2 4 ) 


B„ = o 


1 e 


IT . 

4 1 


2 V 3 ( G - c ) 5/4 


J 


where, in deriving the second row, asymptotic expressions (B.12) and 
(B.19) for $2 and have been used. Relationship (A. 24) is again in 
qualitative agreement with its classical equivalent as it is in the case 
of decaying solution The classical equivalent is (ref. 21) 
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(B . 25 ) 



1 

c T 

r 

$4 (T) 

Cl. 

dT ^ 


- ^ 


2 . . . 3/2 

3 (ia o x) 


— oo —oo 


with the same nomenclatures as before and with the Hankel function of 

( 2 ) 

the second kind H ^ . The asymptotic expression of (A. 25) for t >> 1 is 


[< 


$4 (t) 


Cl. 


-5/4 

t exp 


4 2 . ,3/2 

e 3 (a o T) 


(B. 26) 


in qualitative agreement with (B.24). 

The set of solutions (4> 2 ' $ 3 ' ^4 > thus constructed has noteworthy 
characteristics for the Wronskian that serves to simplify the analyses 
to follow; 


III 


*2 *3 *4 
*2 *3 *4 

$ 2 *3 $ 4 


= 1 


(B. 27) 


This formula is easily checked by noting the following relationship 


_ T 


d 4>. 


dn dn 

holding for l; 0, 1, 2. 


$3 d 4>, 

-4 dn + — 

w dn 


* 6 
/2 , + 12 

«r 2 “ 


Solution 

Since all the homogeneous solutions have been exhausted the fourth 
solution of Eq. (B.2) needs to be sought from particular solutions of 
the equation with = 1; 

- ilu-c)^ + iuO JL - i ( B . 2 8 ) 

The method of variation of constants provides a workable tool also here, 
since we have all the homogeneous solutions 0 ^ through 0^ at hand. 
According to the theorem cited no differential equation need to be solved, 
and a manipulation using key property (A. 27) leads to the following from 
of the solution; 
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= i$2 W I-j dn - 


w i 2 dn + i<J> 4 f 


1= 4 W dn ( i ; 2 , 3 ) 

i \ 


where n = n >> 1 is a point beyond which the asymptotic expressions for 
<J> 2 through i> 4 are valid. Taking derivatives successively we have. 


, U. n 

1 - i 2 f 

I — — 1 TS 


X, n 

d $3 r 

W 1 3 dn - i — — U 1 1 2 dn 

dn } Q 


l n 

d $ 4 f 

+ i — tt- W dn + i 6 . 


(I = 1,2,3) 


a relationship easily checked using (A. 26), (A. 20) and (A. 21). 

The function 4^ in form of (A. 29) is shown to tend to a slowly 
varying function for n 1 under the same conditions as invoked in 
deriving asymptotic expression (B.19) for $ 2 ; the resulting expression 


<I> ~ -(u-c) (u-c) dn 


(n > n >> l) 


For comparison the classical equivalent of this function is noticed; 


= T + O(otR) 


Cl , in 


■ i } ci , out = (u_c) f (u " c)_2dy + 0(ft2) 
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where t has been defined in (B.13). Coincidence of leading terms of our 
asymptotic expression with classical outer solution is observed also here. 


In view of (A. 31) we can easily show the following relationship to 
hold regarding the Wronskian formed by the four solutions through 


IV 


lW III “ 1 


(B . 34 ) 


Higher Order Approximation 

As described in the main text the full solution of the Or r- Sommer f eld 
equation is constructed through a successive approximation starting with 
Eq. (B . 2 ) whose solutions are now at hand. In practice, the n-th order 
correction is obtained by solving inhomogeneous equations (112) , 

following the same procedure as has led Eq. (B.28) to solution (B.29). 

We have, then, 


$ 


(n) 


Wl 3 j;»>dn - 4 3 


wi 2 J i (n) dn 


+ 


wj| n) dn 


(B. 35) 


with 


J (n) E 2$ (n 15 - i 

i i 


(u-c) 1) dn - 


r n 


. (n-2) 

dn 


(B. 36) 


Utilizing, in these formula, the same approximation as employed in 
deriving asymptotic expressions for 0^ and $ 2 , we 9 et 


(n) 


rn 


(u-c) 


dn 


n (u-c) J n 


f (u-c)o| n 1 ^dn. (n > n >> l) (b.37) 

J n 
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his expression , being substituted into (115) with (i;l,2) 

jiven by (B.32) and (B.19), respectively, is readily seen to have exact 
coincidence with Heisenberg’s outer (inviscid) solutions (ref. 20). The 
series solutions for and are provied to be uniformly convergent if 
the lower bound y (= en) is taken such that u-c varies only slowly for 
y < y. In fact, then, the series sum up to yield 

F 1 - -(u-c) sinh A(y-y) (B.38) 

F 2 - (u-c) cosh A(y-y) (B.39) 


Asymptotic Boundary Condition 

It is advisable from convergence relief point of view to replace the 
boundary condition of type (102) at a far boundary (y = y) by an equiv- 
alent asymptotic one; 

AY + Y f « 0 (B . 4 0 ) 

This is permissible if, at y = y properly chosen, 


A 


i« 

jj 

u-c 


> O(aR) 


(B.41) 


a condition accessible to most of the laminar flows and the wall-turbulent 
flow discussed in the main text. Imposing (B.40) instead of (102) enables 
us to choose a smaller value for y in favor of convergence in the 
successive approximation. 


Comparison with Existing Stability Calculations 

Although the method described here is designed primarily for wall- 
turbulent flows, we have invoked no assumptions to rule out its applica- 
bility to the classical (laminar) solutions. In fact, the classical case 
is simply a limiting case (ii 0) in Eq. (106) of the present formalism. 
The whole set of laminar solutions forms, in this sense, a subset of 
solutions achieved by the present method. Reproduction, therefore, of 
some of the classical stability characteristics will serve for checkout 
of the program. 

Figure Al shows the neutral stability curve for Blasius ' flow, where 
the comparison is made with those by Lin (ref. 21) and Wazzan et al. 
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(ref. 29) obtained through the classical asymptotic method and purely 
numerical one, respectively. The classical asymptotic method is subject 

-/3 

to inaccuracy of O(otR) , whereas the present calculation includes error 

4 2 

of 0(u ) which is caused by truncating successive approximation at a . 

These limitations are reflected on the behavior of each curve deviating 

from the computational one with a decrease of R and an increase of a, 

respectively. Agreement, however is considerably improved in the present 

case. Experimental data due to Schubauer and Skramstad (ref. 20) are 

also shown for comparison. Their seeming agreement with Lin's curve at 

lower Reynolds numbers is due to a fortuitous reason: Good agreement is 

not to be expected for lower Reynolds numbers where the stability curve 

becomes sensitive to factors arising from nonparallel streamlines of the 

Blasius flow. It is properly corrected by Saric et al . incorporating 

these elements (ref. 31) . A more precise comparison among the three 

methods are seen in Figure A2 , where calculated longitudinal fluctuation 

| Y ' | is shown together with those by classical asymptotic method due to 

Schlichting (ref. 32) and by the computational method (ref. 29), also 

with the experiment (ref. 30). Prediction by the present method is seen 

to have the best correlation with the measurement in the case tested. 


Euler's Transformation for Improved Convergence 


Our main objective of computing the wall-turbulence stability charac- 
teristics faces a difficulty not encountered in the laminar case; the 
problem of convergence in successive approximation (109). It appears 
that the lower bound rj of the integral in (B.36) was not chosen so as to 
secure uniform convergence. Since, however, the actual computation has 
shown that expansion (109) proves to be a (complex) alternating series, 

one can apply Euler’s transformation to speed up convergence. The trans- 

00 

formation is well known as a means of finding the sum I a of an 

n=0 n 

alternating series which is not absolutely convergent, a tool invented 
by Euler (1755) and applied in its full length to problems of the bound- 
ary layer by Meksyn (ref. 33). Operation of the transformation converts 
the sum to 


00 


l 

n=0 


a 

n 


00 


l 

n=0 


b /2 (n+1) 
n 


with 


n 


= l 


n! 


r ^ Q r ! (n-r) I 
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securing much faster convergence owing to alternating structure of (a r ). 
The method applied to sum an asymptotic series is utilized here: First 

we find a having the minimum absolute value 

a = min . { a } 
m n 

If m = 1, the Euler transformation is applied to the whole series. If 

m > 1, we apply the transformation only to the partial sum; { a^ } with 

n >. m + 1. Repeated use of relevant alternatives at each stage is con- 

, ,N i 

tinued for a given finite sum of N = 10 terms until b x „/Z b settles 

N n 

down within a prescribed error bound. 
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Lin (1945) 



Figure Al. Neutral Stability for the Blasius Flow: Comparison of the Classical Asymptotic Method 

(Lin, ref. 21), the Numerical Method (Wazzan et al., ref. 29) and the Present Method 
wit:: experiment due to Schubauer and Skramstad (ref. 30) . 
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Figure A2 . The Longitudinal Fluctuation, or the First Derivative of the 
Eigenfunct ion Y' Computed or Measured at (R = 2080, a = .307) on the 
Upper Branch of the Blasius Stability Curve: Comparison of the 

Classical Asymptotic Method ( Schlichting , ref. 32) , the 
Numerical Method (Wazzan et al., ref. 29) and the 
Present Method with Experiment due to 
Schubauer and Skramstad {ref. 30) . 
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